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Abstract: Fourier transforming the virtual photon transverse momentum in 
7* (</_!_) + N — > f processes allows new insight into hadron dynamics as a function of 
impact parameter b. I discuss how previous analyses of charge density based on elastic 
and transition form factors (/ = N, N*) can be generalized to any multi-hadron final state 
(/ = irN, irirN, DA C , . . .). The 6-distribution determines the transverse positions of the 
quarks that the photon couples to, and can be studied as a function of multiplicity, the 
relative transverse momenta, quark masses and polarization. The method requires no fac- 
torization nor leading twist approximation. Data with spacelike photon virtualities in the 
range < Q < Qmax provides a resolution Ab > 1/ Q m ax in impact parameter. 
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1. Charge density from elastic form factors 

The quark density of a hadron in transverse (impact parameter) space b is given by a two- 
dimensional Fourier transform of the electromagnetic form factor of the hadron [1, 2, 3, 4]. 
For a proton with helicity A the density distribution is 

Pq/Ar(fe) = / |^e" i ^^ T (P + ,iq ± ,A|j + (0)|P + ,-iq ± ,A) 

= l Q C0 ^QMbQ)F 1 (Q 2 ) (1.1) 

where F\(Q 2 ) is the Dirac form factor. It is not immediately obvious why precisely this 
definition corresponds to a density. Until recently it was in fact common to define the charge 
density in terms of a i/iree-dimensional Fourier transform. Such a definition is, however, 
not compatible with relativistic effects. Quarks in the proton move with nearly the speed 
of light, v q ~ c, so a photon cannot give a sharp picture of the charge distribution at an 
instant of time. However, the transverse velocity = p q /E q decreases with the energy 
Eq of the quark. In the Infinite Momentum Frame (IMF) = and a high resolution 
picture of the charge distribution can be obtained in the transverse plane, i.e., as a function 
of impact parameter as in (1.1). Formally, the IMF is equivalent to quantization at equal 
Light-Front (LF) time, x + = t + z. More intuitively, a photon moving along the negative 
z-axis interacts at fixed x + . The Fourier transform in (1.1) is in fact defined in a frame 
with photon momentum q + = 0. It is also important that the matrix element in (1.1) 
involves only the j + component of the quark current, and that the initial and final states 
have opposite transverse momenta. 

To see why p(b) merits being viewed as a charge density one needs to expand the 
hadron h state in terms of its quark and gluon Fock components taken at equal x + [5], 

n,Aj i=l V * j j 

x Vnfci, k ±i , Xi) \n; XiP + , XiP ± + k ±i , A;) (1.2) 

For a proton the Fock states n would include \uud), \uudg), \ uuduu), . . ., the infinite and 
complete sum of all quark and gluon states, integrated over the longitudinal momentum 
fraction xi and the relative transverse momentum of each parton 1 . The unique property 
of this LF expansion is that the wave functions ip!^(xi, Aj) do not depend on the hadron 
momentum P + , P± . Hence the same wave functions ?/>n describe the initial and final states 
in (1.1). Each parton i carries a share X\ of the parent hadron's longitudinal and transverse 
momentum, and a relative transverse momentum k± i . 

1 This is a formally exact expansion, but possible contributions from partons with Xi = (zero-modes) 
will be neglected. 
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When the initial and final states in (1.1) are expanded in their Fock states (1.2), the 
impact parameter distribution of a quark q is found to be [3] 

n,Xi i=l i i 

x \^( Xi ,bi,Xi)\ 2 Y, e k5 {2) (b-b k ) (1.3) 

k 

where the wave functions ^(xj, 6j, Aj) are related to the momentum space wave functions 
ip^(xi, Aj) of (1.2) by standard (two-dimensional) Fourier transforms of the Thus 
Pq/h(b) indeed is a charge density: the probability that there is a quark k in the hadron at 
impact parameter = b (relative to the parent hadron), weighted by its charge e^. 

The usual parton distributions f q /h( x iQ 2 ) measured in hard inclusive processes may 
likewise be expressed in terms of the LF wave functions, 

rk±<n d 2 k± . 



/qAfoM 2 ) = E[n / dx * I " ^^l^^-E^^E*^) 



n,\i i=l 



x \^(x h k ±i ,\i)\ 2 S{x k 



(1.4) 




(a) (b) 

Figure 1: (a) The elastic proton form factor measured by the ep —> ep scattering amplitude. In 
the frame q + = its Fourier transform (1.1) over gives the charge density in impact parameter, 
(b) The quark distribution in longitudinal momentum f q /h( x ,Q 2 ) (1-4) may be obtained via QCD 
factorization in the Q 2 — > oo limit from Deep Inelastic Scattering, ep — > eX. 

It is apparent that the quark distributions in impact parameter (1.3) (Fig. la) and in lon- 
gitudinal momentum fraction (1.4) (Fig. lb) give similar and complementary information 
on hadron structure. However, there are also important differences. First of all, the ex- 
pression (1.4) is not exact due to the neglect of the "Wilson line", indicated by the vertical 
Coulomb gluon exchanges in Fig. lb. This contribution arises from the rescattering of the 
struck quark in the color field of the hadron and adds coherently to the bound state wave 
function. There is no Wilson line in the form factor since it has only one photon vertex. 
Hence the expression (1.3) of the impact parameter distribution p q /h(b) in terms of the LF 
target wave functions is formally exact. 
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Parton distributions like f q /h( x , n 2 ) are obtained from hard inclusive scattering through 
QCD factorization at leading power ("twist") in the Q 2 — > oo limit, and depend on the 
factorization scale fi. Corrections to the hard subprocess of higher order in a s must be 
taken into account in the extraction of f^/hip -,^ 2 ) from data. Conversely, the expression 
(1.3) for the impact parameter distribution only has corrections from higher orders in the 
QED coupling a. Gluon corrections to the photon vertex are included in the sum over 
Fock states n 2 . Data in the whole range of Q 2 is used, in fact the integral in (1.1) is over 
all spacelike q 2 = — q ± 2 = — Q 2 < 0. A finite range < Q < Q max = q±max limits the 
resolution in impact parameter to Ab > l/Q ma x- 

2. Charge distribution of inelastic processes 

The analysis described above has been applied to data on elastic (7*iV — > N) and tran- 
sition (7*iV — > N*) form factors [6, 7]. The sum over initial and final Fock states in 
(1.3) remains diagonal for transition form factors, but the product of wave functions 
ipn ( x i-> bii ^iT^Pn ( x «> ^ii ^i) i s no longer positive definite. The impact parameter distribu- 
tion nevertheless reflects the transverse positions of the quarks which couple to the virtual 
photon. 

Here I shall describe the generalization of the impact parameter analysis to any pro- 
cess 7*i — > /, where i and / are arbitrary (multi-)hadron states [8]. The possibility to 
study how the impact parameter distribution depends on the type and relative momenta 
of the produced hadrons allows new insight into hadron dynamics. I shall use the process 
r y*N(p) — > ir(pi)N(p2) to illustrate the procedure. 

The Fock expansion (1.2) is valid for any hadronic state, and thus also for \ir(pi)N(p2)) ■ 
However, we need to make sure that irN states with different total momenta pf = p + q = 
Pi +P2 are described by the same LF wave functions. The obvious guess is to parametrize 
the hadron momenta in the same way as the parton momenta in (1.2): 

Pi = xpt p 1± = xp f± + k± 

1 (2.1) 
p+ = (l-x)p+ p 2± = (1 -x)p f± - k ± 

We may specify the trN state using a wave function ^ (x, k±) of our choice, 

\*N(p+,p f± ;9t)) ee ^ ^ J ^^(x,k ± )\n(pi)N( P2 )) (2.2) 

where x and k± are independent of the total irN momentum pf. The asymptotic |7riV) state 
(2.2) then has an LF expansion of standard form, with quark and gluon wave functions 
given by ty* and the Fock state wave functions an d °f (1-2)- It should be noted 
that the wave functions which determine the density distribution are the ones at the time 
of the photon interaction (x + = 0), not those of the asymptotic (x + — > oo) irN state. The 

2 The renormalization of vertex corrections would introduce a scale if the charge density were factorized 
into a photon vertex and wave functions. However, such scale dependence is absent in the measurable 
density (1.3). I am grateful for a discussion with Jian-Wei Qiu on this point. 
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|7rJV) state evolves with x + : the hadrons fly apart at large times and converge toward the 
photon vertex as x+ — > 0. As the pion and the nucleon get close to each other they start 
to interact and can form resonances. Hence the A(^f*N — > irN) amplitude has a dynamical 
phase, unlike the spacelike form factor A(j*N — > N) which is real. A single hadron in 
the final state has a stationary time development, i.e., its Fock state wave functions are 
independent of x+. 

The impact parameter analysis of the 7* A?" — > ttN transition amplitude is similar to 
that of form factors. In the frame where 

p= (p + ,p~,-igj_); q = (0 + ,q~,q ± ) ; p f = (p + ,p~ + q~, \q ± ) (2.3) 

the Fourier transform of the j+ current matrix element can be expressed as a diagonal sum 
over Fock states, 

AfN{b) " / ^-^ b ^f(Pf)\j + mN{ P )) = ^£[n jf **i J we* 

x 5(1 - ^ x t )5 2 £ Xi&^r^, bi) ^2e k 5 2 (b k -b) (2.4) 

i i k 

where iph are the LF wave functions of the state (2.2) at the photon vertex (x+ = 0). 
AfN gives the impact parameter distribution of the quarks to which the photon couples, 
when the center-of-momentum of the initial (TV) and final (/) state is at Xjfoj = 0. 

Since the amplitude {f{pf)\ j + (0)\N(p)) has a dynamical phase due to final-state in- 
teractions the Fourier transform in (2.4) generally requires a partial wave analysis. Alter- 
natively, the Fourier transform of the square of the amplitude, 



S fN (b)^ J ±-(f( Pf )\j+(0)\N(p)) 2 = j d 2 b q A fN {b q )A) N {b q -b) 

(2.5) 

gives a convolution of the impact parameter distributions (2.4). Now b is the difference 
between the impact parameter of the quark struck in the amplitude and in its complex 
conjugate. This also gives a measure of the transverse distribution of the active quarks. 
The quantity Sfwib) has an imaginary part if the squared amplitude in (2.5) is asymmetric 
for — > — The imaginary part thus measures the azimuthal correlation between q± 
and a transverse direction defined, e.g., by a relative transverse momentum between the 
particles in the final state (such as k± in (2.1)). 

The square of the amplitude is obtained from the measured cross section - with the 
caveat that the analysis concerns only the matrix element of the j + component of the 
current. This component dominates at high lepton energies, or may be identified via a 
Rosenbluth separation. 

The above method can be illustrated using the Born level QED amplitude for scattering 
on a muon, j*(q) + fi{p) -> /x(pi) + 7(^2), 

A^' + h„ \- 9 r\ e - ' k± e_ ■ [k± - (1 - x)q ± ] \ 

+ 1+1 [Q±) ~ V 1 (1 - x) 2 m 2 + k ± 2 (1 - x) 2 m 2 + [k ± - (1 - x)gj 2 / 1 ° J 
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where e_ • k ± = e~^ k \k±\/V2. The muons and the photon are taken to have positive 
helicities, the final state momenta are parametrized as in (2.1) and the wave function 
*f(x,k±) of (2.2) is a ^-function in x and k±. The Fourier transform (2.4) gives 



T 6 z (b) = Ki(mb) 



(1 - x) 2 m 2 + k ± 2 2V2it 1 



k± ■ b 

exp l —1 



1 — x 

(2.7) 

In the first term of (2.6) the virtual photon interacts with the initial muon, and this term 
contributes to (2.7) at the impact parameter b = of the target. In the second term 
the virtual photon interacts with the muon after the emission of the real photon, and its 
fo-dependence agrees with the known jjl — >■ jjL + 7 QED wave function. 

If the Fourier transform is taken of the square of the QED amplitude (2.6) as in (2.5) 
the result is 

S^' + hb-x kA-Ae 2 xi k±2/2 5 (2) (b) im \ k ^ cos ^-(pk) Kjjmb) 

\^ {b ' X ' k±) ~ ^ X \[(l-x) 2 m 2 + k ± Y {b) ' m (l-x)W + fc ± 2 2vr(l-x 



Ko(mb) — \mb Ki{mb) ~\ f k± ■ b 



x) 

The three terms within { } correspond, respectively, to the virtual photon interacting (i) 
with the initial muon in both and (A^)* , (ii) once with the intial and once with the 
final muon, and (m) twice with the final muon. The imaginary part can be seen to arise 
from the angular correlation between the lepton scattering plane (defined by b) and the 
relative transverse momentum k 1 in the final state. 



3. Discussion 



The analysis presented here generalizes previous work on charge densities of (transition) 
form factors. Being applicable to any final state in j*N — > f it opens up a new window on 
the dynamics of lepton- nucleon scattering. Data at all spacelike photon virtualities < Q < 
Qmax are used, with an expected impact parameter resolution Ab > l/Q m ax- The absence 
of a QCD factorization removes uncertainties related to the leading twist approximation and 
the factorization scale. The possibilities to study how the impact parameter distribution 
depends on properties of the final state (multiplicity, relative momenta, quark masses, . . . ) 
can give new insight into hadron dynamics. For example, the dimensional scaling observed 
[9] in deuteron photodisintegration, -yd — > pn at 9cm = 90° , suggests transversally compact 
configurations of the deuteron and nucleons. A Fourier transform of the electroproduction 
process, j*d — > pn, can reveal the transverse distribution of the active quarks. 

The absence of QCD factorization also implies less predictions. The Q 2 -dependence of 
the quark distributions / q /Ar(x, Q 2 ) measured in DIS can be calculated, and the universality 
of the distributions tested in other hard processes. The LF wave functions ip n which 
determine the impact parameter distribution in (2.4) are also universal, but are more 
difficult to reconstruct from measured data. The impact parameter analysis discussed 
here thus is complementary to the traditional analyses of hard inclusive (and exclusive) 
scattering processes. 
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